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Abstract 

We examine the cavity resonance tuning of high-Q silicon photonic crystal het- 
erostructures by localized laser-assisted thermal oxidation using a 532 nm continuous 
wave laser focused to a 2.5 fim radius spot-size. The total shift is consistent with the 
parabolic rate law. A tuning range of up to 8.7 nm is achieved with ~ 30 mW laser 
powers. Over this tuning range, the cavity Q decreases from 3.2x10^ to 1.2x10^. Nu- 
merical simulations model the temperature distributions in the silicon photonic crystal 
membrane and the cavity resonance shift from oxidation. 



1 Introduction 

Photonic crystal nanocavities are increasingly employed in photonic studies and applica- 
tions because of their high quality factor (Q) to modal volume (Vm) ratios [HEIIH]. These 
nanocavities are also used in photonic devices of increasing complexity where high accuracy 
of the resonant wavelength is critical. However, due to fabrication imperfections, resonances 
will often deviate from their desired precise values. Several post-fabrication tuning tech- 
niques have been proposed and demonstrated to address this issue. These methods can 
be divided into 2 groups: global and local. Global tuning creates a uniform change over 
the entire chip, [H El IH] whereas local tuning only changes a small area such as a single 
nanocavity [3 13 O [TTl [121 US] • Global tuning is useful for correcting uniform errors but 
cannot address the random local errors that often occur during fabrication. Local tuning 
can be very important to applications such as the all-optical analog to electromagnetically 
induced transparency [2] and optical buffers, which require precisely coupled cavities. An- 
other example is in solid-state cavity quantum electrodynamics where local tuning can be 
used to spectrally match a cavity resonance to a single exciton transition. 




Figure 1: (a) Illustration of laser-assisted local thermal oxidation on a silicon double- 
heterostructure cavity, (b) SEM image of a double-heterostructure cavity. 

In this work, we study the selective thermal oxidation of silicon photonic crystal mem- 
branes with a highly localized laser beam at ambient conditions in order to finely tune the 
resonances of high-Q nanocavities (Fig. 1). Laser-assisted local oxidation is advantageous 
to previously demonstrated tuning techniques because it allows for automation of the tun- 
ing process by using a computer-controlled stage, shutter or optical modulator and in-situ 
monitoring system. Such a system would enable the post-fabrication fine-tuning of large 
numbers of nanocavities. 

Laser-assisted local tuning has been demonstrated with the photodarkening of chalco- 
genide films on GaAs photonic crystals [U] with a Q of ~ 8000. Laser-assisted local tuning 
has also been studied in the oxidation of GaAs photonic crystal L3 cavities [7] with a Q 
of 1800. In this work, we study the precise local tuning in silicon, of photonic crystal 
double-heterostructure nanocavities with high-Qs of ~ 300,000 or higher. 

Continuous-wave lasers have been used in oxidation studies of silicon and silicon-on- 
insulator substrates [HI US], including temperature independent contributions to silicon 
oxidation from photon flux [H]. Using a diffraction limited beam, Deustchmann et al. [T6] 
was able to oxidize lines as narrow as 200 nm at a power of approximately 15 mW. Such 
spatial confinement is possible in thin single crystal silicon films because increased phonon 
scattering reduces heat flow in the lateral direction ^17j . In addition to cavity tuning, local 
oxidation might be applicable to other post-fabrication tasks such as tuning the dispersion 
[THl dU] and surface states ^011^. 

2 Local oxidation cavity resonance tuning 

The photonic crystal double-heterostructure nanocavities p| used in this work were fab- 
ricated by high quality photolithography and dry etching on silicon-on-insulator (250 nm 
thick) substrates, with 117 nm hole radii and 410 nm lattice parameter. In the cavity region, 
the lattice parameter increases to 415 and 420 nm. The waveguide-to-cavity separation is 
6 layers of air-holes. Approximately 1.5 ^m of oxide beneath the photonic crystal region 
was removed as described in Ref. 4. 

A 532 nm diode-pumped solid-state laser, with a collimated beam and power controlled 
using a variable neutral density filter to ~ 20 mW powers, was used. A 60x objective lens 
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(NA of 0.65) focuses the laser onto the chip and is also used for imaging. The spot-size 
is measured using the knife-edge technique. The full-width half-maximum beam waist was 
2.5 /im, corresponding to a maximum energy density at the chip surface of 1x10^ W/cm^. 
The oxidation is carried out in ambient conditions of 22 °C and 20% relative humidity. 

Cavity resonance transmission and radiation measurements were performed as described 
in Refs. 2 and 4. Oxidation results in a blueshift in the cavity resonance as shown in Fig. 
[2] (a) . This is the result of a larger decrease in refractive index from the silicon consumed 
outweighing a smaller increase in refractive index from the oxide generated. Following the 
initial anomalous oxidation, subsequent tuning follows a parabolic rate law where the cavity 
resonance blueshift is observed to be proportional to the square root of the oxidation time 
(Fig. [2] (b) and (c)). The decreasing oxidation rate is attributed to the longer diffusion 
time of oxygen through the thicker oxide in order to reach the silicon-oxide interface where 
the oxide growth occurs |22] . Corresponding radiation measurements of the high-Q cavity 
are shown in Fig. [3| 
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Figure 2: (a) Experimental results showing the blueshift from tuning using an initial laser 
power of 19.5 mW (at the device surface). The same device is then further tuned at 22.5 
mW. The inset shows measurements of the loaded quality factor as the cavity is tuned, (b) 
Fitting of the resonant wavelength shift to the square root of oxidation time for incident 
power of 19.5 mW and (c) for 22.5 mW. 

Oxidation is expected to occur on both the top and bottom surfaces of the silicon mem- 
brane as well as the surface of the air holes because both the 1 /e absorption depth of the 
green laser (1 fim) and the depth of focus of the beam after the objective lens (400 nm for 
60 X objective) are larger than the 250 nm thickness of the silicon membrane. The oxidation 
profile remains highly localized with AFM studies showing lateral dimensions 40% smaller 
than the spot diameter [15j. This can be attributed to a number of factors. The ther- 
mal conductivity of silicon is dominated by phonon transport with a smaller contribution 
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Figure 3: Radiation measurements at three different tuning increments corresponding to 
the numbered positions in Fig. 2. Left inset: near-infrared radiation pattern for the Fabry- 
Perot modes on the left. The dotted hnes indicate the position of cavity waveguide (upper) 
and input /output waveguide (lower). Right inset: near-infrared radiation pattern for the 
cavity mode on the right-side of the spectra. 

from free charge carriers. While the thermal conductivity of silicon is relatively high at 
room temperature, it decreases significantly at higher temperatures with increased phonon 
scattering [23', '24]. In addition, as the substrate thickness decreases from bulk dimensions, 
phonon-boundary scattering increases and the presence of air holes will further decrease 
the phonon mean free path |17j . 

3 Transient effects from oxide surface chemistry 

During the oxidation process, as the cavity is heated to high temperatures by optical absorp- 
tion of the focused green laser beam, there will be a large redshift in the cavity resonance. 
The redshift is attributed to the thermal-optic effect in silicon causing a temperature- 
dependent refractive index change. When the green laser is turned off, the cavity temper- 
ature quickly returns to room temperature and the large redshift disappears. But within 
1-2 minutes after oxidation (i.e. after the green laser is turned off), a smaller magnitude 
redshifting of the cavity resonance is observed. The rate of this redshifting is rapid at first 
and decreases over time, taking many hours to reach a stable value. This transitory effect 
is caused by water molecules on the surface of the cavity. 

The surface chemistry effects occur both during and after laser irradiation, and are 
illustrated in Fig. [4| Upon heating the cavity there is a temporary blueshift from oxide 
surface and oxide bulk dehydration. After the cavity cools back down, there is a slow 
redshift caused by gradual rehydration of water molecules onto the cavity. The total shift 
ranges from ~ 100 pm to hundreds of picometers and appears to be dependent on the 
oxide thickness. Laser exposure at low powers (less than 1 mW for 10 minutes) shows a 
completely reversible blueshift, indicating no real oxidation has occurred (i.e. there is only 
a temporary blueshift from cavity dehydration followed by a rehydration redshift returning 
the cavity resonance to its original wavelength). 

Silicon dioxide terminated by hydroxyl groups (SiOH) is hydrophillic and will readily 
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Figure 4: Transitory surface chemistry effects, (a) Water molecules absorb onto the oxide 
surface at room temperature and desorb from the oxide surface at elevated temperatures, 
(b) During laser irradiation, the cavity experiences a resonance blueshift from oxide dehy- 
dration in addition to silicon oxidation. After the cavity cools, water will slowly rehydrate 
the oxide surface resulting in a gradual redshift. 



adsorb water molecules, as schematically shown in Fig. |4] (a). At temperatures above 
^ 170 °C the hydrogen bonded water molecules will desorb. The process is reversible but 
at higher temperatures (400 °C) the hydrogen in the hydroxyl groups can sometimes be 
removed resulting in a hydrophobic siloxane surface. The long times required to complete 
the redshift are indicative of a slower diffusion- limited process [25\ . While water absorp- 
tion is known to be much less significant for thermally grown oxides than deposited oxides 
[27] , studies of laser grown oxides indicate their composition are less dense than traditional 
thermally oxidized films due to presence of suboxides, especially for thinner films [28]. In 
order to obtain reliable results, our measurements were taken either immediately after ox- 
idation (less than 1 min) when the cavity region was still dehydrated or immediately after 
the cavity was re-heated at sub-oxidation threshold powers. 



4 Thermal oxidation 

Prior to oxidation, the test chip had a native oxide (approximately 12 A in thickness) as it 
was exposed to ambient conditions. In addition to the native oxide, we observed that initial 
oxidation can occur at very low laser powers. Because of this, an accurate initial oxidation 
threshold power was difficult to determine. This initial growth was anomalous in respect 
to the Deal-Grove and Massoud models for oxide growth [22^ I29j . It has been shown by 
x-ray photoelectron spectroscopy studies that the initial ~ 22 A of oxide growth occurs at 
a very high rate [_30j. At this stage, even at very low powers oxidation might occur. 

Following the initial anomalous oxidation, subsequent tuning follows a parabolic rate 
law. The oxidation resonance tuning is permanent and stable. Repeated measurements 
over the span of many days shows no change within the measurement error. The uncer- 
tainty in the cavity resonance measurements is approximately it 100 pm and primarily 
attributed to random thermal fluctuations in the ambient environment. This uncertainty 



5 



is the limiting factor in determining the minimum tuning increment. A number of other 
factors will affect the minimum achievable tuning increment including: power stability, ex- 
posure time accuracy, beam targeting accuracy and beam shape uniformity. Another group 
has demonstrated a resonance tuning increment of 100 pm using atomic force microscope 
induced nano-oxidation of GaAs 




Figure 5: (a) SEM image of cavity after local oxidation tuning at laser power of 35 mW. 
The white ring indicates oxide charging effects during SEM imaging. The center region is 
darker because of slight melting. Inset: Hole melted through the silicon membrane after 
irradiation at ~ 35 mW for several minutes (different device from the main figure), (b) Local 
oxidation tuning over a larger wavelength range. The blue region (upper-left) corresponds 
to the data shown in Fig. 2 (a). The upper-right arrow corresponds to the SEM image in 
Fig. 1 (b). The inset shows measurements of the loaded quality factor. 

The inset to Fig. [2] (a) shows quality factor measurements during the oxidation process. 
During the initial oxidation there is a sizeable drop in the Q. Afterwards the Q appears to 
remain stable upon further oxidation. The initial drop in Q might be attributed to increased 
optical scattering from the interfacial layer between the silicon and oxide. Surface states 
might also result in increased absorption at the interface. In this work, the local oxide was 
left on the chip after the tuning. Removing the oxide with an HF dip [7J has been shown 
to improve the Q and will also further blueshift the resonance. 

The decrease in overall air-hole size can be used to estimate the oxide thickness. The 
SEM image in Fig. [T] (b) shows the double-heterostructure cavity after local oxidation 
resulting in a 2.1 nm cavity resonance blueshift. Image analysis of the experimental device 
indicates a hole size change smaller than the error of the SEM measurements (it 5% of the 
hole diameter). In addition, a redshift of 0.8 nm was observed in the cavity resonance after 
SEM imaging the cavity region (Fig. [5](b)). This is likely caused by surface contamination 
of the sample in the SEM chamber during imaging. It is well known that organic molecules 
from vacuum pump oil can be deposited onto the substrate surface by the focused electron 
beam during SEM imaging [8]. For this reason, SEM imaging was generally avoided while 
the sample was being tuned. Fig. [s] (a) shows an SEM image of the sample after cavity 
resonance tuning over the full range shown in Fig. ^ (b). The total resonance blueshift is 
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8.7 nm. Note that if the SEM-induced contamination redshift did not occur, then the total 
resonance blueshift would be 9.5 nm. Using hole-size analysis it is estimated that 1 nm of 
resonance blueshift corresponds to 1.9 nm of oxide growth. A comparison with calculated 
resonance shifts is provided in the next section. 

The measured loaded Q of the cavity over the entire tuning range decreases from 3.2 x 10^ 
to 1.2x10^ as shown in the inset of Fig. [H] (b). It has been shown by 3D FDTD analysis 
that, for the double-heterostructure cavity design used, an increase in hole radii of ~ 12% 
corresponds to a decrease in Q of ~ 80% [Tj. In our experimental case, our hole radii 
increases by ~ 8% and Q decreases by ~ 63% (only considering the silicon; note that the 
overall hole radii actually decreases because of the oxide growth). Hence, the sensitivity of 
Q to the hole size plays a significant role in the observed drop in Q during tuning. 

5 Numerical analysis 

Three-dimensional finite-element simulations (COMSOL Multiphysics) are used to estimate 
the temperature increases caused by optical absorption at different laser powers. The model 
uses the static heat equation: 



VKr)Vr] = -[l-i?(T)]/oa(T)exp -^-^ exp - / a{T{x,y, /))dz' ] (1) 



where k{T) is the temperature-dependent thermal conductivity, R{T) is the reflection co- 
efficient and a{T) is the temperature-dependent absorption coefficient. The incident CW 
laser beam has a Gaussian spatial distribution with a 1/e-squared beam spot-size a and a 
center intensity Iq. The reflection coefficient takes the form R{T) = Rq + c{T — Tq) [STl [T1] . 
The temperature dependent absorption is experimentally fitted as a{T) = aoexp(r/ro) 
[32j. Both the absorption and thermal conductivity of silicon are temperature dependent. 
The absorption of silicon increases significantly at higher temperatures. The thermal con- 
ductivity for bulk silicon at room-temperature is 148 WK~^m~^ but as the temperature 
increases, the thermal conductivity will decrease. For a thin porous photonic crystal slab, 
the thermal conductivity can be reduced even further. In-plane thermal conductivity is re- 
duced to approximately 90 WK~^m~^ at room temperature for silicon-on-insulator (SOI) 
devices with a silicon thickness of 260 nm |33], mainly due to increased phonon boundary 
scattering [24^ l3^ l35] . Periodic structures on single-crystalline silicon membranes can have 
thermal conductivity values as small as ~ 6.8 WK^^m"^ due to coherent phononic effects 



In order to determine the thermal conductivity for our device, we used our finite-element 
model along with experimental measurements. Silicon melting is observed to occur at 
incident laser powers of ~ 35 mW as shown in Fig. [5](a). This was used as a reference point 
for our model and the resulting room-temperature thermal conductivity was estimated to be 
67 WK~^m~^. The simulated temperature profile around the double-heterostructure cavity 
is shown in Fig. [6] (a). Cross-sections of the temperature profile corresponding to incident 
power of 30 mW and 35 mW are shown in Fig. [6](b) along with the laser intensity profiles at 
those powers. The temperature profiles show how the temperature gradient increases due to 
higher absorbed energies and lower thermal conductivity at higher membrane temperatures. 




[36] . 
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The local maximum temperature for different incident laser powers is shown in Fig. [6] (c). 
The features observed in Fig. [5] (a) are consistent with results of the numerical model Fig. 
[6] (a). 

(a) (b) (c) 




Figure 6: (a) Finite-element simulation (COMSOL Multiphysics) of temperature distribu- 
tion across silicon double-heterostructure cavity during 532 nm laser irradiation at 35 mW. 
(b) Solid lines represent the temperature distribution as a function of distance from the 
center of the laser beam. Dotted lines represent the intensity profile of the laser beam, (c) 
Simulation results of local maximum temperature versus laser power. Temperatures range 
from room temperature to the melting point of silicon. 

The finite-difference-time-domain (FDTD) method [37] with sub-pixel averaging [38] is 
used to calculate the effects of silicon oxidation on the cavity resonance (Fig. 7). The 
native oxide is assumed to be 1.2 nm. We also assume that if 1 nm of silicon is oxidized 
(consumed), the resulting oxide grown is 2.17 nm (0.46 ratio). As the silicon is consumed, 
the total thickness of the slab (silicon -|- oxide) will increase, while the air-hole radii will 
decrease. Simulation results in Fig. [T] (b) show that every 0.6 nm of oxide grown results in 
1 nm of resonance blueshift. By comparison, the previously discussed SEM hole analysis 
estimates that 1.9 nm of oxide growth results in 1 nm of resonance blueshift. There are 
several possible explanations for this discrepancy. It is possible that the amount of oxidation 
on the bottom-side of the silicon membrane is less than the top-side if the beam focus is less 
than optimal. In addition, the laser grown oxide might be less dense than furnace grown 
thermal oxide [2^. Finally, the slight melting of the silicon membrane at 35 mW laser 
irradiation might have resulted in a decrease in the air-hole size without a proportional 
increase in oxide growth. 

6 Conclusion 

We have demonstrated the tuning of high-Q double- heterostructured silicon photonic crystal 
nanocavities using laser-assisted local thermal oxidation. Cavity Q decreases from 3.2x10^ 
to 1.2x10^ over the range of oxidation times and laser powers examined. The effects of 
water absorption and thin oxide growth were also observed. Numerical simulations were 
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Figure 7: (a) Calculated electric field Ey profile of high-Q mode supported by double- 
heterostructure cavity, (b) Calculated wavelength shift of the resonant mode due to local 
oxidation of the silicon photonic crystal membrane. 



used to model the temperature distribution in the silicon photonic crystal membrane and 
resonance shift of the optical mode due to oxidation. 
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